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We investigate theoretically the dynamics of squeezed state generation in the nonlinear systems 
possessing the transition from regular to chaotic dynamics in the limit of a large number of photons. 
As an example, the model of kicked Kerr oscillator is considered. We show the direct correlation of 
the degree of squeezing and the value of local Lyapunov instability rate in corresponding classical 
chaotic system. 
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In recent years the problem of generation of the squeezed states of electromagnetic field has attracted a great 
attention of both the theoreticians and experimentalists fl]-^| . As a rule, in the experiments on squeezing the quantum 
fluctuations are small compared to the mean value of the field modes involved into a nonlinear interaction . Already 
for a rather long time it was mentioned that a light squeezing can be increased near the bifurcation points between 
different dynamical regimes The reason for such enhancing consists in a strong diverging of one of the field 

quadrature component and corresponding squeezing of another one at the vicinity of instability threshold |^J^] . Among 
with well studied parametric media , the increase of squeezing was predicted for the interaction of the field with 
2-level atoms inside a high-Q cavity at the dynamical regime corresponding to the separatrix ||||. However, all 
papers devoted to the study of the enhanced squeezing due to the appearance of dynamical instabilities in optical 
\^ ■ systems dealt with the integrable or near integrable systems with regular dynamics [To[ . 
£SJ \ In this paper, we consider the squeezed light generation by quantum nonintegrable optical system obeying the 
transition from regular to chaotic dynamics in the classical limit. As for other problems of quantum chaos fll2] , |l3|] , 
we consider the semiclassical limit when a great number of quantum levels N 1 are involved into the dynamics. 
Our prime goal is to demonstrate the sufficient increase of the degree of light squeezing at transition to quantum 
chaos. We also found direct correlation between the degree of local instability in corresponding classical system and 
the degree of light squeezing. 

In spite of our consideration is valid for any nondissipative quantum system with one and half degrees of freedom, 
we will demonstrate our main results on the enhanced squeezing for some particular model of the nonintegrable optical 
system - the nonlinear oscillator periodically forced by the classical field. In the interaction picture the Hamiltonian 
^ ■ has the form (h = 1) 

A 

H = Ab+b + -b +2 b 2 + eN^ 2 (b + b+)F(t), (1) 

: 2 

where Bose operators b and b + ([b,b + ] = 1) describe a single-mode of the quantum field and A is proportional to 
third-order nonlinear susceptibility of a nonlinear medium. The last term in (Q) corresponds to a coupling of the 
oscillator with external classical periodically modulated field containing a large number of photons N 1 (e is a 
coupling constant, F(t) is a periodic function of time and A is a detuning of the mode frequency from the carrier 
frequency of the external field). 

The Hamiltonian ([j]) describes, for example, a high-Q cavity filled by a medium with Kerr nonlinearity and excited 
by an external laser field fli}] . This effective Hamiltonian may also govern the interaction of a laser field with a 
high-density exciton system in a semiconductor 15]. The diffcren variants of discussed model are very popular in 
both quantum optics and quantum chaos |l3|. 17-pfj|| studies. 

To investigate the nonlinear dynamics of the model and the dynamics of squeezing, it is natural to use the cumulant 
expansion technique (2lj as a variant of the general 1/iV-expansion method p2]| . Recently it was realized that both 
methods could be applied not only to the integrable [^|2^j2^] but also to the nonintegrable quantum systems |l8| , |2"3|J 
In this paper, we use well-adopted for the problems of quantum optics th 1 /iV-expansion method suggested in 

Introduce new normalized operators for annihilation and creation of photons as a — b/N 1 / 2 , a+ = b+fN 1 / 2 with 
commutation relation [a, a + ] = l/N. In the classical limit (N — > oo), one has two commuting c-numbers. Now the 
Hamiltonian (|l) can be rewritten as H — NH^, where has the same form as (^) with an account of replacement 
b — > a, b^ — ► Qj and A — * g = XN. It may be shown that the dependence g ~ N correctly performs the time scale of 
energy oscillation for Kerr nonlinearity in the classical limit [E5[ . 
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Let initially the field is in the coherent state | a) = exp{Naa + — Na*a) | 0) corresponding to the mean photon 
number ~ N. From the Heisenberg equations for a, a 2 and their Hermitian conjugated equations, and adopting the 
normal operator ordering, we have the following equations of motion for averages over coherent states 

i^-(a) = {V), i^-{{8a) 2 ) = 2{VSa) + (W), i^{8a*8a) = -{V*Sa) + (Sa*V), (2) 

lib (JLL (JjV 

where V = |ff , W = jf-g^r, {{8a) 2 ) = {a 2 ) - (a) 2 , {8a* S a) = (a+a) - (a+)(a). However, the set of equations (§) is 
not closed and actually is equivalent to the infinite dynamical hierarchy system for moments and cumulants ]7|pl|] . To 
truncate it up to the cumulants of the second order, we made the substitution a — > {a) + 8a, where at least initially 
the mean (a) ~ 1 and the quantum correction | 8a(t = 0) |~ TV -1 / 2 <C 1. Using the Taylor expansion of the functions 
V and W, we have, in the first order of 1/N and after some algebra, the following self-consistent system of equations 
for the mean value and the second order cumulants 

iz={V) z +q, (3a) 
+ ,m B, (3b) 



da J „ V da 



where B = {Sa*8a) + -jy , C = {{8a) 2 ), z = {a) and subscript z means that the values of V and its derivatives are 
calculated at mean value z. Involved to the equation ( ^a| ) the small quantum correction q ~ 1/N has the form of the 
second differential of V as follows 

l a „ , 1 fd 2 V\ „ l(d 2 V\ „ / d 2 V \ / I 



* - i dV l*= 2 {^) z C+ 2{^) z C * + {da^) z [ B 2N 

The initial conditions for the system (^) are -B(O) = (7(0) = and some arbitrary z(0) which is of the order of 
unity. 

It is easy to see that the equations of motion ( pb| ) and ( pc[ ) for the cumulants can be obtained from the classical 
equations by linearization near z (substitution z — ► z + A«, | 5a |<c| 3 |), if one writes the dynamical equations for the 
variables (Aa) 2 and | Aa | 2 . But there still exists the principal difference between the linearization of the classical 
motion equations and the equations for quantum cumulants (y|): it is impossible to obtain the initial conditions for C 
and B from only initial conditions for linearized classical equations. 

Define the general field quadrature as Xg = aexp(—i8) + a + exp(i6>), where 9 is a local oscillator phase. A state is 
said to be squeezed if there is some angle 9 for which the variance of Xg is smaller then the variance for a coherent state 
or the vacuum 0,^]. Minimizing the variance of Xg over 9, one can determine the minimum half-axis of the quantum 
noise ellipse. Then, the condition of the principal squeezing is S = 1 + 2N{{\ 8a | 2 }— | {{8a) 2 ) |) = 2N(B— | C \) < 1 
. The determination of the principal squeezing S is very useful because it gives the maximal squeezing measurable 
by the homodyne detection O. 

Let us now compare the dynamics of the principal squeezing for regular and chaotic motion. Initially Gaussian wave 
packet spreads when it propagates through nonlinear medium. But it still exists the time interval of the well-defined 
quantum-classical correspondence during which the wave packet center follows path in the phase space governed by 
semiclassical equations of motion (|3a|). Moreover, because our equations of motion for cumulants (|3b| ) and (^), in 
fact, coincide with the equations arising in the definition of the maximal Lyapunov exponent |]l2| , |l3|| , we can apply 
simple physical arguments on strong deformation of the classical phase volume at chaos for prediction of the strong 
squeezing of the noise ellipse at quantum chaos in the semiclassical limit. For chaotic dynamics, the distance in the 
phase space between two initially very closed trajectories D grows exponentially with time D{i) ~ exp(At), where 
A is the maximal Lyapunov exponent. Due to the presence of local strong (exponential) local instability inherent 
for underlying classical chaotic dynamics, a quantum noise ellipse may be strongly stretched in one direction and 
squeezed in another direction. As a result, the value of the principal squeezing in average exponentially decreases in 
time at chaotic dynamics. The stretching and squeezing of noise ellipse at quantum chaos is much stronger than for 
the case of regular and stable dynamics, when the distance between two initially closed trajectories in phase space 
increases in time power-wise resulting in only power-wise decreasing of the principle squeezing in time. 



2 



In order to illustrate this general picture, we return to Hamiltonian ([!]) and choose the form of F(t) as a periodic 
sequence of kicks: F(t) = J2^ = _ oc S(t — nT), where S(t) is the Dirac 5-function. In an experiment, a sequence of short 
light pulses can be generated by a mode locked laser. Now, by means of standard technique Q], we obtain from the 
differential equations (0) the coupled maps for the mean value and for the second order cumulants. In this paper, we 
present only the form of the resulting map for the mean value p7|j2^ | 

zn+i =e- iT ( A+ ^-^ 2 )(z n -ie) + q n (C n ,B n ,z n ), (4) 

where the subscript n indicates the value of a function just before the action of the n-th kick. The use of the maps 
instead of the differential equations sufficiently simplifies the computation and reduces numerical errors. 

In the classical limit N — > oo and q — > 0, the map (Q) possesses the transition from regular to chaotic dynamics. 
Fig. 1 illustrates the behavior of the classical map: the phase portrait and the time-dependence of intensity z n 2 
for regular (Fig. la,b) and chaotic (Fig. lc,d) dynamics. 

Define the distance in the phase space between two initially very closed trajectories as D{t) — [(5x) 2 + (Sy) 2 ] 1 ^ 2 , 
where x = Rez and y = Imz. The Fig. 2 shows the time-dependence of the logarithms of the principal squeezing 
S (Fig. 2a), the normalized distance between two initially very closed trajectories d — ND (Fig. 2b), and the 
normalized value of quantum correction Q = Nq (Fig. 2c) for large but finite number of quanta N = 10 9 . In this 
figure, curve 1 corresponds to the regular motion, and curves 2 and 3 - to chaotic motion with slightly different values 
of the Lyapunov exponent. As it is evident from this figure, the most strong local instability determines the highest 
degree of squeezing. It should be noticed that the difference in the magnitude of principal squeezing for chaotic and 
regular motion may achieve of several orders during only several kicks. 

As follows from ([|), for chaos, both cumulants B and C increase also exponentially in average, resulting in an 
exponential growth of the quantum correction q. All our analysis is valid, if and only if the values of the quantum 
correction and of the second order cumulants are much less than the mean values: q, |C|,.B<^|.z|~l. Under the 
conditions of chaos, this means that the equations (|^) are correct only during the time interval t <f ~ A~ x In N. The 
time interval t* coincides with time scale determining the well-defined quantum-classical correspondence in chaotic 
systems jL7[. Nowadays the time scale t* , during which classical chaos is revealed as a quantum transient, is established 
for the different models of quantum chaos |p^| , ^8|j24| . Moreover, only during this time interval it is possible to define 
the Lyapunov exponent for a quantum system |29| . Our numerical experiments on the model of kicked Kerr oscillator 
demonstrate, that for the chaotic motion at gT = 7, A/g — 1 and e = 0.1, the value of the principal squeezing 
S for N = 10 7 is practically indistinguishable from its value for N — oo up to the six kicks. The increase of the 
number of photons N results in the corresponding increase of the time interval for applicability of our description of 
squeezed dynamics. On other hand, the number of photons N > 10 7 initially pumped to the system is quite realistic 
for contemporary experiments on light squeezing Q . 

In summary, we have showed that among with very narrow class of the 

integrable systems near threshold of an instability, there is another wide class of potentially effective systems for 
the enhanced light squeezing - the systems with quantum chaos operating during the time scale of the well-defined 
quantum-classical correspondence. 
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FIG. 1. The nonlinear dynamics of kicked classical nonlinear oscillator: phase portrait and time-dependence of the intensity 
\ z„ | 2 for the regular behavior at gT = 3 (a, b) and for the chaotic motion at gT = 10 (c, d). Initial condition is zq — 1 and 
e = 0.1, A/s = 1. 
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FIG. 2. The time-dependence of the principal squeezing (a), of the distance between two initially closed trajectories (b), and 
of the quantum correction (c). Curve 1 corresponds to the regular dynamics at gT = 3, curve 2 - to the mild chaos (gT = 7) 
and curve 3 - to the hard chaos (gT — 10). The average photon number is N = 10 9 , initial condition and other parameters are 
the same as in Fig. 1. 
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